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ABSTRACT 

Galaxy-galaxy mergers and close interactions have long been regarded as a viable mechanism for channeling 
gas toward the central supermassive black holes (SMBHs) of galaxies which are triggered as active galactic 
nuclei (AGNs). AGN pairs, in which the central SMBHs of a galaxy merger are both active, are expected 
to be common from such events. We conduct a systematic study of 1286 AGN pairs at z — 0.1 with line-of- 
sight velocity offsets Av < 600 km s"' and projected separations < 100 /i^q kpc, selected from the Seventh 
Data Release of the Sloan Digital Sky Survey (SDSS). This AGN pair sample was drawn from 138,070 AGNs 
optically identified based on diagnostic emission line ratios and/or line widths. The fraction of AGN pairs with 
5 h^Q kpc ^ rp < 100 /z^o kpc among all spectroscopically selected AGNs at 0.02 < z < 0.16 is 3.6% after 
correcting for SDSS spectroscopic incompleteness; ~'30% of these pairs show morphological tidal features in 
their SDSS images, and the fraction becomes > 80% for pairs with the brightest nuclei. Our sample increases 
the number of known AGN pairs on these scales by more than an order of magnitude. We study their AGN and 
host-galaxy star formation properties in a companion paper. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei - quasars: 
general 



1. INTRODUCTION 

Binary supermassive black holes (SMBHs) are thought to 
be a generic outcome of galaxy-galax y mergers in the hierar- 
chical paradigm of structure formation (iBegelman et al.ll980l : 
iMilosavljevia & Merrittl200TtlYul20 02'). give n that most mas- 
sive g alaxies harbor central SMBHs ( Kormen dy & Richstond 
ll995URichstone et alJ[T998l) . If both SMBHs accrete mate- 
rial during the same stage of a galaxy-galaxy merger, they 
can be identified as a pair of active galactic nuclei (AGNs). 
However, it is difficult to predict from first principles when 
one and in par ticular both SMBHs become active during 
a merger (e.g., Shlosman et al."199^ lArmitage & NatarajanI 
^002; Wada 2004; Dotti et al. 2007j). The frequency of AGN 
pairs can constrain models invol ving galaxy m erger rates and 
tidally-triggered AGN activity (lYu et al.ll20Tli) . Their host- 
galaxy properties may offer clues to the external and internal 
conditions under which both SMBHs can be activated during 
galaxy-galaxy tidal encounters. 

When two galaxies merge, the two SMBHs in-spiral, form 
a bound binary, and harden, before they finally coalesce 
dBegelman et al.|[l980h . The coalescence of binary SMBHs 
is expected to be the strongest source of gravitational waves 
in the universe, the detection of which would directly test gen- 
eral relativity a nd could probe SMBH populations out to the 
early universe ( Thomel 19871: lHaehneltlll994l iHolz & HughesI 
l2005h . However, there is still Uttle direct observational ev- 
idence for small-separation (< 10 pc) gravitationally bound 
binary SMBH s on Keplerian orbits (see e.g., the review of 
IColpi & Dottil2009l) . 0402H-379 is the only unambiguous case 
known. This object was detected by Very Long Baseline In- 
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terferometry as a pair of flat-spectru m radio point sources 
with a proje cted separation of —7 pc dRodriguez et al.ll2006l: 
IValtonenet al. 2008). At the early stages of a merger, large- 
separation (kpc to tens of kpc) AGN pairs where the galac- 
tic potential dominates can readily be spatially resolved at 
cosmological distances. They provide boundary/initial condi- 
tions for small-separation pairs. By characterizing the large- 
separation population, we may gain indirect but useful infor- 
mation on the small-separation population. 

There are an increasing number of examples of interacting 
AGN pairs with kpc to tens-of-kpc sep arations. Such objects 
have been i dentified in the X-ray (e.g.,'Komossa et al."2003; 
Ballo et alJl2 004; Bianchi et al. 2008; Piconcelli et al. 201(1 
Mazzarella et al. 2011) the radio (e.g., Owen et al. 1985), or 
the optical (e.g., B arth e t al. 2008; Comerford et al. 2009b] 
iGreen et"ani2010t [Liu et a l. 2010a). However, most of the 
known cases of kpc-scale AGN pairs have been serendipitous 
discoveries. To systematically identify and characterize the 
population of AGN pairs, we selected 167 obscured AGNs 
with double-peaked narrow emis sion lin es from a parent sam- 
ple of 14,756 AGNs (Liu et al J 1201 Obi) in the Sloan Digital 
Sky Survey (SDSS; York et all f 
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lution near-IR imaging and spatially resolved optical spec- 
troscopy, we identified four kpc- scale A GN pairs amo ng 43 
of the 167 double-peaked AGNs dLiu et al. 2010a; Shen et alj 
1201 lb . While our work has demonstrated that the double-peak 
approach allows us to systematically identify kpc-scale AGN 
pairs, it has two inherent limitations. First, the approach is bi- 
ased against pairs with line-of-sight (LOS) velocities smaller 
than ~'150 km s"' due to the limited spectral resolution of 
SDSS spectra. Second, it is biased against pairs with sepa- 
rations > a few kpc because pairs with projected separations 
larger than 3" (the diameter of the SDSS fiber) will not fall 
within a single SDSS fiber 

To circumvent the limitations of the double-peak method 
and to build up the sample of AGN pairs along the merger se- 
quence, we take a complementary approach to identify them. 
We start with a sample of spectroscopically identified AGNs 
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and select physical pairs with special attention to those that 
are in mergers or close interactions. Because each AGN in 
a pair has a separate SDSS spectrum, the current approach 
is not limited by the SDSS spectral resolution or by the 3" 
fiber size. The current sample is sensitive to an earlier stage, 
with separations of a few kpc to tens of kpc, of galaxy-galaxy 
mergers than is the double-peak approach. Such systems are 
quite scarce, so a statistically meaningful study requires ho- 
mogeneous and large-area redshift surveys such as the SDSS. 
While galaxy pairs with projected angular separations smaller 
than 55" cannot both be observed spectroscopically on the 
same SDSS plate due to the finite diameter of the fibers, those 
in the overlap region s of adjacent plates can both be spectro- 
scopically observed j Strauss et al.l l2002: Blanton et al. 2003). 

The present paper is the first in the series in which we 
present a statistical sample of AGN pairs from the Seventh 
Data Release (DR7: lAbazajian et alJi200^ of the SDSS. In 
^we describe sample selection, and discuss incompleteness 
and selection biases. We address the frequency of AGN pairs 
among the parent sample of optically identified AGNs in |3] 
We discuss the implications of our results in Ml and conclude 
in ^ In a companion paper dLiu et al.ll201 laT hereafter Pa- 
per II), we will characterize host-galaxy recent star forma- 
tion and BH accretion properties of the AGN pair sample 
presented here. Throughout we assume a ACDM cosmology 
with n„ = 0.3, r^A = 0.7, and Hq = 70 hjo km s"' Mpc"'. 

2. SAMPLE SELECTION 

In this section we describe our approach to select AGN 
pairs. In 32.11 we construct a parent AGN sample optically 
identified from the spectroscopic sample of SDSS DR7. We 
then draw AGN pairs with LOS velocity offsets Av < 600 km 
s"' and projected separations < 100 /i^q kpc ( ^2.2\ . To mit- 
igate contamination due to pairs that are closely separated but 
are not interacting and to focus on AGN pairs that are unam- 
biguously experiencing strong tidal encounters, we visually 
inspect the SDSS images of all AGN pairs and further iden- 
tify a subset, which we call the "tidal" sam ple, that exhibits 
unambiguous morphological tidal features ( 32.3t . We discuss 
our sample completeness and selection biases in 32.41 

2.1. The Parent AGN Sample 

Our parent AGN sample includes objects drawn from the 
SDSS "naiTow-line "^ AGN sample (Ka uffmann et al. 2003a^, 
broad-line AGN s (iHao et al.1 l2005blal). narrow-hne quasars 
dZakamska et al.' '2003^, 'Reyes et alj |2008|), and broad-line 
quasars (Schneider et al. 2010), all from the DR7^. Below we 
describe the selection criteria for each population. The red- 
shift cut for all objects is 0.02 < z < 0.33. The lower redshift 
limit is to balance the need to include nearby galaxies and to 
avoid redshifts where peculiar velocities can cause substantial 
deviations from Hubble flow; the upper redshift limit is to en- 
sure Ha coverage in the SDSS spectra for AGN identification. 

1. Narrow-line AGNs. First, we select AGNs from 
the MPA-JHU SDSS DR7 galaxy catalog^ (see 

' The "narrow-line" classification here is not meant to be st rict. Some 
AGN s in this sample show broad-line components as well (e.g.. IHao et alj 
l2005bl) . 

* The narrow-line quasar sample of IReves et all J2008ft was selected from 
a sample corresponding to ~80 % of the D ata Release 6 (DR6) spectroscopic 
database ( Adelman-McCarthy e t al.l200^ so we may miss a small number of 
narrow-line quasars post DR6. 

' http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 1. — Redshift distributions of the parent AGN sample (black solid 
curve) and the subsamples it contains (color curves). We use "NLAGN" to 
denote narrow-line AGNs according to the Kauffmann et al. I2003a |) crite- 
rion, "BLAGN" for broad-Une AGNs according to the Hao et al. (20053) 
classification that are not contained in the MPA-JHU DR7 AGN catalog, 
"NLQSO" for narrow-line quasars in the Re ves et al.i (i2008i) sample but not 
in the MPA-JHU DR 7 AGN catalog, "QSO" for quasars in the catalog of 
Schneider et al." EOHJ), and "Faint QSO" for quasars above the flux limit of 
the SDSS main gala xy sampl e but whose luminosities are too low to be in- 
cluded in the iSchneider et alj i2010i) catalog. Refer to ^for details. 

lAihara et al.l 1201 U for a description of the cat- 
alog) drawn from the SDSS DR7 spectroscopic 
database; this database includes 869,580 targets'^ spec- 
trally classified as galaxie s by t h e specBS p i peline 
dAdelman-McCardiy et al.1 l2008t lAaaraetal.1 I2OIII; 
see also discussion in Blanton et al. 2005) or quasars 
that are target ed as galaxies (Strauss et al. 2002] 
lEisenstein et akl l200lh and have redshifts z < 0.7. 
We adopt redshifts z and stellar velocity dispersions 
cr* from the specBS pipeline. Additional spectro- 
scopic and photometric properties such as emission- 
line fluxes and stellar masses are taken from the 
MPA-JHU data product. The emission line measure- 
ments a re from Gaussian fits to continuum-subtracte d 
spectra (iBrinchmann et al.ll2004t iTremonti et al.ll2004b . 
The MPA-JHU emission-line fluxes have been nor- 
malized to SDSS r-band photometric fiber magnitude 
and have been co rrected for Galac tic foreground ex- 
tmction follow ing lO'Donnelll ([1994) using the map of 
ISchleeel et al] (119981) . The stellar mass estimates are 
total stellar mas ses derived from population synthe- 
sis fits using the iBruzual & C hariot (2003 ) models to 
SDSS broadband photometry flKauffmann et al.ll2003bl; 
[Salimet al. 2007). 

We have two selection criteria: (1) emission lines H/3, 
[Oiii]A5007, Ha, and [Nii]A6584 must afl be de- 
tected with signal-to-noise ratios (S/N) > 3, and (2) the 
diagnostic emission-line ratios [O iii] A5007/H/? and 
[Nii]A6584/Ha suggest that the dominant excitation 
mechanism i s from AGN rathe r than stellar photo- 
ionization (Baldwin et al. 1981; Osterbrock & Poggel 
119851; IWUeux & Osterbrock ■1987b according to the 
empirical criterion, 

^"g^^Q"^^/"^^> log([NIlVHa)-0.05 "^-^' 

as suggested by 'Kauffm ann et al] ( 12003 a*). We address 
the effect of the adopted S/N threshold for emission- 
line detections on our result s in ^ 129,277 objects 
meet these two criteria. The iKauffmann et alJ (l2003ah 

Here and throughout, we quote the numbers of unique galaxies after 
excluding duplicate observations of the same galaxy. 
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empirical criterion includes AGN-H n "composite" ob- 
jects, which represent over half the sample. If we were 
to adopt inste ad the theoretic al "starburst limit" sug- 
gested by Ke wlev et alj <2001h . 



log([OIII]/H/3)> 



0.61 



log([NII]/Ha)-0.47 



+ 1.19, 



(2) 



in order to exclude the AGN-H n composites, the 
sample would be reduced to 50,624 objects. We in- 
clude composites in our analysis, because the results 
would be otherwise biased against starburst-dominated 
systems with moderate AGN activity. In a study 
of hard X-ray selected moderate-luminosity AGNs at 
z < 0.05 from the Swift Bu rs t Alert Tel e scope AGN 
sample (iTueller et al.l 1201 Ol) . iKoss et all (1201 Oh find 
that 33% of AGNs in merging systems show com- 
posite or star-forming diagnostic lines. Similarly 
iGoulding & Alexanded (2009) found evidence of star 
formation in the optical emission lines of AGNs iden- 
tified based on Spitzer infrared spectra of IR luminous 
galaxies. 

2. Narrow-line quasars. Second, we supplement the par- 
ent sample with 1 12 objects from the narrow-line 
quasar sample of iReves et al.l (120081) . These objects 
satisfy our selection criteria for narrow-line AGNs but 
are not already contained in the MPA-JHU SDSS DR7 
galaxy catalog because they were not part of the main 
spectroscopic galaxy sample. 

3. Broad-line AGNs. Third, we supp lement the sample 
with broad-line AGNs selected by iHao et al.1 (l2005bh 
extended to SDSS DR7. These are obj ects selected 
from the SDSS main galaxy sample dStrauss et al.l 
|2002) at z < 0.33 which have a broad Ha component 
(FWHM > 1200 km s"') from multi-Gaussian fits with 
a rest-frame Ha equivalent width (EW) > 3 A. There 
are 9573 such objects at 0.02 < z < 0.33, and 3710 of 
them are included in the MPA-JHU DR7 AGN cata- 
log or in t he S DSS DR7 quasar catalog (see below). 
iHao et all (l2005b) show that the narro w-line compo- 
nents o f most broad-line AGNs satisfy the lKewlev et alj 
(1200 ll) criteria for narrow-line AGNs. Given the pur- 
pose of this work, we supplement the parent AGNs 
with the new 5863 broad-line AGNs, and keep the 3710 
overlaps in our "narrow-line" AGN category without la- 
beling them more carefully. As shown by Reyes et al.l 
(l2008i) . it can be unclear whether an AGN is a narrow- 
line or a broad-line object when the broad component 
is relatively faint or when the forbidden emission lines 
have non-Gaussian profiles (such as extended wings 
and/or double peaks). The classification can be fur- 
ther complicated by the presence of scattere d broad-line 
region light in obscu red quasars (e.g., Zakamska et al.l 
l200aiLiultani2009l) . 

4. Broad-line quasars. We next supplement the parent 
AGNs with 2944 objects at 02 < z < 0.33 from 
the SDSS DR7 quasar catalog (ISchneid.er et al .1120101) 
The SDSS quasar catalog of ISchneideret'aiT(l2010l) 
has a luminosity cut M, < -22.0. There are quasars 
wh ose luminosities are to o low to be included in 
the ISchneider et al.l (1201 Oh quasar catalog but have 
extinction-corrected Petrosian 14.5 < r < 17.77 mag 



3.0 



2.5 



m 

E 2.0 



1.0 



0.5 



I I I I I I I I I I I I I I I I I I I I I 



IS.". 



I I I I I I I I I I 



lull 



20 



40 



(h; 



60 80 

; kpc) 



100 200 400 
§ of AGNs 



Fig. 2. — Projected separation rp vs. LOS velocity offset Ai' (on a loga- 
rithmic scale) for AGN pairs (black dots and dotted histogram) and the sub- 
set with identified morphological tidal features (blue open circles and solid 
histogram). The Ai' distributions have been corrected for spectroscopic in- 
completeness ( j|2.2.U . We caution that values below Av = 30 km s"' are 
dominated by redshift measurement uncertainties. 

consiste nt with the flux lim it of the SDSS main galaxy 
sample dStrauss et al.ll2002h . To include those as well 
we searched through all objects targeted as "quasars" 
(iRichards et al.ll2002h at 0.02 < z < 0.33 in the spectro- 
scopic sample of SDSS DR7 and further supplemented 
the sample with 173 broad-line objects after rejecting 
stellar contaminants. These quasars are not included in 
the DR7 version of the broad-li ne AGN sample as th ey 
were not targeted as "galaxies" (IRichards et al.ll2002l) . 

The resulting parent AGN sample consists of 138,070 op- 
ticafly selected AGNs at 0.02 < z < 0.33. We Hst the num- 
bers of galaxies in each selection category in Table |2] Among 
the parent AGN sample, 94% are included in the MPA-JHU 
DR7 galaxy catalog, 2. 9% of which are c lassified as broad- 
line AGNs according to lHao et"aLl (l2005bh because they have 
weak FWHM > 1200 km s"' Ha components. This catalog 
serves as the starting point from which we draw the "pair" 
sample based on projected and radial separation criteria, and 
a subset which shows tidal features in SDSS images. To put 
our results on AGN pairs into context and to compare with 
"ordinary" AGNs, we also draw control AGNs from this par- 
ent sample (Paper II). Figure[T]shows the redshift distribution 
of the parent AGN sample and its sub-populations. 

2.2. AGN Pairs 

We select AGN pairs from our parent sample of AGNs 
based on the following criteria: (1) projected separation rp < 
100 Z;,,) kpc and (2) LOS velocity difference Av < 600 km 
s"'. We choose 100 kpc to minimize contamination due to 
pairs that are closely separated but are not interacting while 
maintaining a large dynamical range in separation. 100 kpc 
is somewhat larger than the t ypical thresh old va lues adopted 
in gal axy pair studies ( e.g., IBarton et al . 2000; E llison et alj 
,2008;iDarg et al.ll2010h . The velocity-offset threshold also 
balances contamination and statistics; it is comparable to the 
galaxy pairwise velocity dispers ion for p rojected separations 
0.15 < rp < 5 /Z7(') Mpc (e.g., IZehavi et al. 2002). While 
galaxy pairs with small projected separations and small veloc- 
ity offsets are likely to be bound and merge within a Hubble 
time, the maximum velocity offset of an interacting galaxy 
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Fig. 3. — SDSS gri-co\oT composite images l ILupton et alj|2004l) of SIX 
examples of AGN multiples. The top row displays three systems of AGN 
triples, whereas the bottom row presents two quadruples and a quintuple. Red 
circles indicate the locations of AGNs. North is upand east is to the left. The 
FOV of each panel is 100" X 100". See Table [T] for coordinates, redshift, 
and other photometric and spectroscopic measurements of each AGN. 
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Fig. 4. — Redshift vs. projected separation (r,,) for objects in AGN pairs. 
The blue dashed curve marks the separation of 55", below w hich th e ob- 
served pair sample is incomplete due to SDSS fiber collision ( j|2.2.1) . The 
black solid histograms show the observed redshift and r,, distributions of 
AGN pairs, whereas the red dashed histograms represent the distribution after 
correction for SDSS spectroscopic incompleteness. The dotted histogram is 
the redshift distribution of the parent AGN sample scaled down by an arbi- 
trary factor for display purposes. The magenta solid curve denotes the sep- 
aration of 2", below which the observed pair sample is incomplete because 
smaller separation pairs are not resolved by SDSS imaging (lAbazaiian et alj 
l2009h . 

pair depends on local environment. To check whether we 
missed a significant population of interacting galaxies with 
even larger velocity offsets, we show the distribution of Av 
of the AGN pairs in Figure |2] The distribution function drops 
sharply at Av > 300 km s"' and the threshold value 600 km 
s"' is on the tail of the distribution, verifying that the Av 
threshold is adequate. 

We visually inspected the SDSS images of all spectroscopi- 
cally selected pairs and rejected 34 false positives due to mul- 
tiple spectroscopic observations of different locations of the 
same galaxy. The final AGN pair sample includes 1286 sys- 
tems with 1244 doubles, 39 triples, two quadruples, and one 
quintuple totalling 2616 AGNs. Figure|3]shows the SDSS im- 
ages of six examples of AGN multiples, although they all have 
large separations without obvious tidal features. Using spa- 
tially resolved spectroscopic follow-up of AGN pairs in the 
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Fig. 5. — Diagnostic line ratios for AGNs in the pair sample (cyan dots) 
and those with tidal features (open circle). Gray scales indicate numbe r densi- 
ties of 31,179 emission-line galaxies from the SDSS DR4 I Kauffm ann et all 
12003d) . The solid curve displays the empiiical separation between H ii re- 
gions and AGNs (Kauffmann et al. 2003a), the dashed curve denotes the the- 
oretical "starburst limit" from Kewlev et al. (2001), and the dotted line repre- 
sents the empirical division between Seyferts and LINERs (Ho et al. 1997|). 
Pure star-forming ("SF") galaxies lie below the solid curve, AGN-dominated 
objects (Seyferts above and LINERs below the dotted line) lie above the 
dashed curve, and AGN-H ii composites ("C") lie in between. 

sample with < lOJz^kpc and having a third close com- 
panion, iLiu et al.l (1201 Ibh discover a kpc-scale triple AGN. 
In Table [l] we list basic photometric and spectroscopic mea- 
surements for all 1286 systems of AGN pairs or multiples. 
Figure |4] shows their distributions of redshift and projected 
separations. While the parent sample spans a redshift range 
of 0.02 < z < 0.33, only 24 pairs are found at z > 0.16 mainly 
reflecting the sparseness of the parent AGN sample. To miti- 
gate uncertainties due to small number statistics at the higher 
redshift end, we restrict our statistical analysis to the 2568 
AGNs atz < 0.16. 

We show in Figure |5] the optical diagnostic line ratios 
[Oiii] A5007/H/3 and [Nii] A6584/Ha of AGNs in the pair 
sample that are contained in the MPA-JHU SDSS DR7 galaxy 
catalog. Sixty-three percent AGNs in the pair sample li e 
above the theoretical "starburst limit" of Kewley et al. (200 1*). 
In Paper II, we will use this sample to examine star formation 
and accretion properties in the host galaxies of AGN pairs and 
their correlations with pair separation. 

2.2.1. Correction for SDSS Spectroscopic Incompleteness 

Our sample of AGN pairs is incomplete for pairs with 
projected angular separations smaller than 55" because both 
galaxies are required to have SDSS spectroscopic observa- 
tions. Due to the finite size of the SDSS fibers, galaxy pairs 
with separations smaller than 55" will not both be observed 
unless they fall in the overlap regions on adjacent plates 
(IStrauss et al.ll2002t iBlanton et al.ll2003h . We correct for the 
fiber collision incompleteness by supplementing the observed 
AGN pairs with (C- l)A^p systems randomly drawn from those 
with separation smaller than 55", where Np is the observed 
number of pairs. Neglecting the tiny fraction of triples and 
higher multiples, the ratio between the intrinsic and the ob- 
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TABLE 1 
SDSS AGN Pairs 
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(mag) 
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f^tidal 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


J00:02:49.07+00:45:04.8 


. . 0388 


345 


51793 


0.0868 


5.9 


9.5 


63 


16.08 
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J00:02:49.44+00:45:06.7 


. . 0685 


593 


52203 


0.0865 


5.9 


9.5 


63 


16.43 


1 


3 


J00:02:57.21+00:07:50.5 


. . 0685 


531 


52203 


0.0901 


45.4 


76.0 


109 


16.15 


2 


2 


J00:02:58.59+00:08:31.0 


. . 0387 


072 


51791 


0.0897 


45.4 


76.0 


109 


16.67 


2 


2 


J00:03:23.74+01:05:47.3 


. . 1490 


321 


52994 


0.0993 


12.2 


22.4 


22 


17.98 


2 


2 


J00:03:23.74+01:05:59.5 


. . 0387 


620 


51791 


0.0994 


12.2 


22.4 


22 


17.42 


2 


2 


J00:04:25.78-09:58:54.4 


. . 0650 


598 


52143 


0.1122 


24.1 


49.1 


53 


16.41 





2 


J00:04:26.66-09:58:34.4 


. . 0651 


392 


52141 


0.1120 


24.1 


49.1 


53 


16.57 


2 


2 


J00:04:31.92-01:14:11.7 


. . 0388 


282 


51793 


0.0887 


42.3 


70.1 


45 


17.22 








J00:04:33.25-01:13:34.4 


. . 1490 


241 


52994 


0.0889 


42.3 


70.1 


45 


18.47 


2 






Note. — The full table is available in the electronic version of the paper Col.(l): SDSS names with J2000 coordinates given 
in the form of "hhmmss.ss+ddmmss.s"; Col. (2): Spectroscopic plate number; Col. (3): Fiber ID; Col. (4): modified Julian date; 
Col. (6): angular separation; Col. (7): transverse proper separation; Col. (8): Hne-of-sight velocity offset; Col. (9): SDSS r-band 
model magnitude corrected for Galactic extinction; Col. (10): flag for sub-populations in the parent AGN sample (refer to ^for 
details). "0" through "2" are for naiTow-line AGNs contained in the MPA-JHU DR7 catalog, where "0" stands for Seyferts, "1" 
for LINERs, and "2" for composites, respectively, according to the Kewley et al. (2001) and Kauffmann et al. i 2003a) criteria for 
separating AGNs and composites from H ii regions, and the Ho et al. 1 1997) criterion for separating Seyferts from LINERs (Figure 
|5). "3" for narrow-line quasars in the Reyes et al. (2008) sample but not in the MPA-JHU DR7 catalog, "4" for broad-line quasars, 
a nd "5 " for broad-line AGNs according to the .Hao et al., (2005b^ selection; Col. (II): flag from visual identification of tidal features 
( j)2.3t . "0" for "non-interaction", "I" for "ambiguous", "2" for "tidal", and "3" for "dumbbell" systems. 




Fig. 6. — SDSS gfi'-color composite images of 36 examples of AGN pairs with tidal features. North is up and east is to the left. The FOV of each panel is 50" 
X 50", except for those in the last row whose FOV is 100" X 100". Objects are arranged in increasing order of transverse proper separation rp from left to right 
and from top to bottom. Binary AGNs in the first, second, third, fourth, fifth, and last row have Vp <5,5 <rp < 10, 10 < rp < 15, 15 < < 20, 20 <rp < 30, 
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Fig. 7. — Examples of ambiguous systems in the pair sample. Objects are 
arranged in increasing order of transverse proper separation r,, from left to 
light and from top to bottom. The FOV of each image tile is 50" X 50" . We 
mark the targets with red boxes when additional objects are within the FOV 
shown. 



served pair fraction is 

( 1 + 2/')/lap/spec,lap + ( 1 + 1 " /lap)/spec 



C(/lap,P) = 



(l+2p)/lap/,pec,lap 



(3) 



where /bp ^ 30% is t he plate overlap fraction on the sky 
jStoughton et alj l2002h . p is the ratio of pairs to isolated 
galaxies, /spec.iap ~ 0.99 is the SDSS spectroscopic com- 
pleteness in the plate overlap region, and /spec ~ 0.92 is 
the spectro scopic completeness in the un-overlapped region 
(iBlanton et al. 2003). C(/iap,p) « l//iap = 3.3 when < 1. 
We show in Figure|4]the distributions of redshift and projected 
separations of AGN pairs after correction for SDSS spectro- 
scopic incompleteness. The corrected redshift distribution of 
the pair sample is similar to that of the parent AGNs, except 
that the pair sample is incomplete for projected angular sep- 
aration < l."4 due to the imaging resolution limit of SDSS 
(the median PSF FWHM of SDSS DR7 images is l."4 in r; 
lAbazajian et al.i 2009). The corrected rp distribution function 
increases with decreasing rp (Figure |4]l. 

2.3. AGN Pairs with Tidal Features 

To mitigate contamination due to pairs that are closely sep- 
arated but are not interacting and to focus on AGN pairs that 
are unambiguously experiencing strong tidal encounters, we 
further select from the AGN pair sample what we dub the 
"tidal" AGN pair sample by visually examining their SDSS 
images for optical tidal features such as bridges, tails, shells, 
or rings. We examine both the gri-coloi composite images 
dLuptonet al.ll2004l) and the cahbrated FITS images in all 
five SDSS bands. We use visual identification as it should 
be more sensitive to low-surface-brightness (LSB) features 
and less subject to false positives than automated techniques 
like model fitting. While quantitative measures of mergers are 
more objective and yield reproducible results (e.g.. lConselicd 
|2b03; Lotz et al. 2004), they may introduce biases toward cer- 
tain merger stage s (e.g., the first p ericenter passage and the 
final coalescence: lLotz et al.ll2008l) . 

We assign each pair a flag of "tidal", "ambiguous", "dumb- 
bell" (discussed below), or "non-interaction" as listed in Ta- 
ble [T] We summarize the numbers of AGNs in each category 
in Table |2] To mitigate uncertainty of the visual inspection, 
one of us (X.L.) classified each object in the pair sample five 
times, with the order of objects scrambled each time and not 



Fig. 8. — Same as Figure|2] but for example dumbbell systems in the pair 
sample. For the systems shown, rp ranges from 2.3 to 9.5 /ly^ kpc. The long 
spike in the upper left panel is due to diffraction from a nearby star 

knowing the results from previous identifications. The identi- 
fication among five trials is identical for 80% of the objects; 
for objects with different classifications from each round (us- 
ing "0" for non-interaction, "1" for ambiguous, and "2" for 
"tidal"; there are few "dumbbell" systems and their classifi- 
cation was almost always the same), we take the median as 
the final result. By inspecting the SDSS single-scan images 
we found 245 AGN pairs with unambiguous tidal features. 
In Figure |6] we show SDSS gri-coloi composite images of 
36 such examples chosen to span the whole range of separa- 
tions. The system shown in the third column and second row 
of Figure |6] is Mrk 266, which is a luminous infrared galaxy 
at z = 0.028. Its northern nucleus is optically classified as a 
Seyfert 2, whereas the southern nucleus is a LINER. Its dou- 
ble AGN nature is confirmed by Chandra X-ray observat ions 
(Mazzarella et al. 2011; see also iBrassington et al.ll200'7l and 
,Wang & GaQ.201Q) . 

2.3.1. "Ambiguous" and "Dumbbell" Systems 

We excluded from our tidal sample those objects we call 
"ambiguous" pairs of AGNs in which both members exhibit 
ordered bar or spiral features but show no clear signs of 
galaxy-galaxy interactions. While these bar or spiral features 
could arise from galaxy-galaxy interactions, they may also 
be induced by internal instabilities. There are 169 "ambigu- 
ous" pairs and we show six examples in Figure]?] We rejected 
these pairs because they are not as convincing as pairs with 
clear merger-induced tidal features (such as disruptive asym- 
metries or tidal bridges). 

In a ddition, we excluded ten "dumbbell" systems (e.g., 
Valen tijn & Casertanolll988h . which have double nuclei in a 
single host galaxy, but show no optical tidal features indica- 
tive of ongoing interactions between the two nuclei. Figure 
[8] shows six such examples. Three of the ten systems contain 
disk components, whereas the others have optical morpholo- 
gies reminiscent of 3C75. 3C75 is a double radio source at 
the center of the galaxy cluster Abell 400 whic h shows evi- 
dence for interactions only in its twi n radio jets (lOwen et al.l 
[1981 11^^1988"; ' Beers et al.|[T992l) . We did not find evi- 
dence for inte ractio ns in the radio by examining the FIRST 
(iBecker et al.H 19951: 1 White et al.lll997l) images of the "dumb- 
bell" pairs, although neither the spatial resolution nor sensi- 
tivity is high enough for a strong constraint. 

2.3.2. Using the Co-added SDSS Stripe 82 Images to Identify 
Low-surface-brightness Features 
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TABLE 2 
Sample Summary 





Selection 


; 

Sub-sample 


Parent 


Pair 


Tidal 


Ambiguous 


Uumbbell 




(1) 




(3) 


(4) 


(5) 


(D) 


/"7\ 


0.02 < z 


< 33 


NLAGN 


129277 


2529 


501 


327 


20 






NLAGN-Kewley 


50624 


1583 


312 


190 


18 






BLAGN 


5564 


79 


9 


10 









NLQSO 


112 


1 


1 












QSO 


3117 


7 


1 


2 





0.02 < z 


< 0.16, observed. . . 


NLAGN 


112242 


2488 


482 


322 


20 






NLAGN-Kewley 


44168 


1555 


299 


186 


18 






BLAGN 


4595 


78 


8 


10 









NLQSO 


25 


















QSO 


279 


2 





1 





0.02 < z 


< 0.16, conected. . . 


NLAGN 


112242 


8116 


2487 


1377 


118 






NLAGN-Kewley 


44168 


5035 


1552 


802 


106 






BLAGN 


4595 


236 


42 


33 









NLQSO 


25 


















QSO 


279 


11 
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Note. — Col.(l): selection criteria. The top block is the basic selection, the middle block is restricted 
to z < 0.16 to mitigate uncertainties due to small number stat istics at higher redshifts, and the bottom block 
is after connection for incompleteness due to fiber collisions ( j)2.2.U ; Col. (2): subsamples in our selection of 
the parent AGN sample. "NLAGN": naiTow-line AGNs according to the Kauffmann et al. r2003a) criterion 
which is our default, "NLAGN-Kewley": narrow-line AGNs according to the Kewley et al. (2001) criterion, 
"BLAGN": broad-line AGNs, "NLQSO": narrow-hne quasars in the Reyes et al. (2008) sample but not in 
the MPA-JHU DR7 catalog, "QSO": broad-line quasars. Refer to 0for details of each p opul ation; Col. (3): 
number of AGNs in the parent AGN sample; Col. (4): number of AGNs in the pair sample ( j|2.2) ; Cols.(5)-(7): 
number of AGNs in the tidal, ambiguous, and dumb bell c ategories according to our visual classification of the 
SDSS images. For details of each category, refer to 32.31 

below the se nsitivity of individual SDSS scan images (e.g., 
ISmirnovaet al. 2010), we examine the co-added images of 
a subset of AGNs that are located in the SDSS Stripe 82^. 
Stripe 82 is a sky region covering ^ 270 deg^ along the Ce- 
lestial Equator in the Southern Galactic Cap, which the SDSS 
has imaged ^100 times, producing co-added imag es ~2 mag 
more sensitive than individual scan images (Aba zaiian et al.l 
l2009h . 

We find that ^ 27% of AGN pairs (6 of 22 pairs) in the "am- 
biguous" category show tidal features in the co-added Stripe 
82 images; the LSB missing fraction is lower 8.5%, 5 of 
58 pairs) for the "non-interaction" category. Figure 0shows 
two such examples. We correct for the LSB missing fraction 
using these estimates of the fraction of the "ambiguous" and 
"non-interaction" cases in what follows. The LSB correction 
factor depends on redshift, since the LSB feature detection is 
redshift dependent due to the limited image resolution. Here 
what we have estimated using the co-added images is the ef- 
fective LSB correction factor averaged over the redshift range 
spanned by the sample, which should be a reasonable approxi- 
mation as the redshift distribution of the Stripe 82 sub-sample 
is similar to the parent pair sample. We supplement the tidal 
AGN pair sample with the 1 1 pairs identified using the co- 
added Stripe 82 images. Our final AGN pair sample with 
tidal features consists of 256 AGN pairs. Table |2] lists the 
numbers of AGNs in each category. The fraction of AGN-H n 
composites with tidal features (out of all composites in AGN 
pairs) is consistent with that of the other AGNs (i.e., narrow- 
line AGNs that satisfy the Kewley et al. 2001 criterion plus 
those in the other AGN categories as discussed in © with 
tidal features (out of all the other AGNs in pairs). Among the 
512 AGNs in the tidal sample, 501 AGNs are contained in the 
MPA-JHU SDSS DR7 galaxy catalog, whose diagnostic fine 
ratios are shown in Figure |5] the other 1 1 objects are selected 
as either narrow-line or broad-line quasars. In the Appendix, 

' http://www.sdss.org/legacy/stripe82.html 
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Fig. 9. — Identifying low surface-brightness tidal features using the co- 
added Stripe 82 images. We show two examples of interacting AGN pairs 
which show tidal features in the co-added Stripe 82 images but not in individ- 
ual SDSS scan images. North is up and east is to the left. The left columns 
are SDSS gr(-color composite images from individual scans, whereas the 
right columns show the co-added r-band images in Stripe 82. The FOV of 
each panel is 200" X 200". The FOV is centered on one of the AGNs in the 
pair, which is connected to the other galaxy via a tidal stream visible in the 
co-added images. 

Our identification for tidal features is limited by the sur- 
face brightness sensitivity and resolving power of the SDSS 
images. The limiting magnitude of individual SDSS scan im- 
ages with a signal-to-noise rat io (S/N) of 5 is r = 23.1 mag 
(in the AB system) for stars (iGunn et al.lll998l) . To quan- 
tify the number of interacting AGNs with LSB tidal features 
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Fig. 10. — Redshift vs. projected separation (rp) for AGN pairs witli tidal 
features. The meanings of the lines are the same as in Figure |4] 
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Fig. 1 1 . — Cumulative fraction of AGN pairs with projected separations 
smaller than r,, (and > 5 kpc). Open circles denote AGN pairs, whereas filled 
circles represent the subset with tidal features. The quoted uncertainties are 
Poisson eiTors. 

we show the SDSS images and fiber spectra for three exam- 
ples of AGN pairs with tidal features. 

2.4. Selection Incompleteness 

The primary factors that limit the sample completeness in- 
clude the resolution and surface-brightness limits of SDSS 
photometry, fiber collisions, the limitations of visual identi- 
fication of tidal features, and the heterogeneity of the parent 
sample. 

Our approach is insensitive to AGN pairs in advanced 
mergers with nuclear separations that are too small to be 
resolved by the deb lending algorithm of SDSS photometry 
jLupton et alj|200l1) . Pairs with nuclear separations < l."4 
will not be included. At z = 0.16, 1" corresponds to 2.7 kpc 
in the assumed cosmology and therefore the sample is incom- 
plete for rp <5 /lyQ kpc. At the low-redshift end, the separa- 
tion limit is set intrinsically by galaxy size, which is a function 



of stellar mass. AGNs in the pair sample have a median stellar 
mass log(M^,/Mp) = 10. 8 using the stellar mass estimates of 
iKauffmann et alj (l2003bh . The z-band half light radius is ob- 
served to be ^^3.3 kpc for this stellar mass (IKauffmann et al.l 
120033) . 

Our tidal sample does not include mergers at very early 
stages before the first peric enter passage; such objec ts have 
no detectable tidal features dToomre & Toomrell 1972b . In ad- 
dition, the strength of the tidally induced features does not 
increase monotonically with decreasing separation when the 
merger advances to later stages; it also depends on the progen- 
itor host-galaxy properties and orbital parameters. We keep 
these separation-dependent limitations and biases in mind 
when discussing our results, although we do not attempt to 
correct for them. Additional populations of mergers with- 
out discernable optical tidal features may inclu de: minor 
mergers whose mass ratios are larger than ~30 dLotz et al.l 
1201 Ol). sphe roidal mergers with large bulge-to-disk ratios 
dMihos & H ernquist 199 ^, and m ergers on significantly ret- 
rograde orbits (Toom re & ToomrellT972D . Some of these bi- 
ases are mass or morphology-dependent and therefore can be 
corrected for using control samples matched in galaxy stellar 
mass (Paper II). 



'0.1 

The fraction of AGN pairs with 5 /i^,' kpc < rp < 100 
/17Q kpc among the parent AGN sample at 0.02 < z < 0.16 
is /p.obs (1265 systems out of 117,141 AGNs) be- 

fore correction for SDSS spectroscopic incompleteness and 
is /p w3.6% after correction. As the pair sample may con- 
tain false positives due to closely separated pairs that are not 
(yet) interacting, this serves as an upper limit for the intrinsic 
fraction of interacting AGN pairs. 

Figure [TO] shows redshift versus transverse proper separa- 
tion and the distributions of these two quantities for AGN 
pairs with tidal features. There are only nine pairs at z > 0. 16. 
This apparent deficiency likely results mainly from the com- 
bination of fewer objects in the parent sample and the de- 
creasing apparent size of tidal features with redshift. The 
redshift distribution of the tidal sample is roughly consistent 
with the parent sample both before and after correction for 
fiber collisions, except that the tidal sample is incomplete for 
smaller separation pairs (projected angular separation < 2"). 
The number of tidal pairs found increases with decreasing rp 
(except for the bin at the smallest scales, which is affected 
by the resolution limit of SDSS photometry), a trend that is 
partly caused by the fact that closer pairs tend to show more 
prominent tidal features. 

We estimate the fraction of AGN pairs with tidal features 
among the parent AGNs at 0.02 < z < 0. 16 as 

(N, + CaNa + C„N„) 

fl = C (/lap , p)/t,obs = C (/lap , p) — , (4) 

where C(/iap,/?) ~ 3.3 is the correction factor for spectro- 
scopic incompleteness (Equation O, A^, is the number of 
galaxies in AGN pairs with tidal features identified in indi- 
vidual SDSS images, Na and A',, are the numbers of galaxies 
classified as "ambiguous" and "non-interaction", Ca = 27% 
and C„ = 8.5% are the correction factors for LSB features 
estimated based on our Stripe 82 experiment (§ 12.3.2) . and 
is the number of AGNs in the parent sample. We find 
/t,ob,s «0.3% before correction for spectroscopic incomplete- 
ness and / wl.0% after correction. Because the tidal sample 
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Fig. 12. — Left: signal-to-noise ratio (S/N) distributions of diagnostic emission-line detection of AGNs in the parent, pair, and tidal samples (after correction 
for spectroscopic incompleteness). Vertical lines indicate the median values for each distribution. Right: dependence of the fraction of AGN pairs on the S/N 
threshold adopted for emission-line detection (for all four emission lines for both AGNs in a pair). Open circles denote AGN pairs, whereas filled circles represent 
the subset with tidal features. The quoted uncertainties are Poisson eiTors. The blue solid curve displays the total number of AGNs in the pai'ent sample as a 
function of S/N threshold. 



is incomplete for interacting AGN pairs without detectable 
tidal features, this serves as a lower limit for the intrinsic frac- 
tion of interacting AGN pairs. 

Figure [TT] presents how our results vary as a function of 
upper limit of r^. Observations of statistical samples of in- 
active galaxy pairs suggest that ^ 30 h^Q kpc is the physical 
scale in projected separation below which galaxy pairs ex- 
hibit sig nificantly hig her star f ormation rates than field galax- 
ies (e . g.,lBarton et al . 2000; L ambas et al.ll2003l:lAionso et alJ 
l2004tlNikolic et al.ll2004) . The fraction of AGN pairs with 5 



hjQ kpc ^ < 30 hjQ kpc is 1 .3%, among which ^ 60% show 
tidal features. 

The fraction of AGN pairs we measure depends on the sen- 
sitivity of our AGN identification. We adopted an S/N of 3 as 
our baseline value for the detection threshold of the diagnos- 
tic emission lines, H/3, [Oiii]A5007, Ha, and [Nil] A6584, 
in the selection of narrow-line AGNs. We now address how 
the S/N threshold affects our results for the fraction of AGN 
pairs. The sample statistics are dominated by narrow-line 
AGNs (Table |2]i, and we thus neglect the effects of changing 
the selection criteria for the other AGN categories. In the left 
panel of Figure [12] we show the S/N distribution of the four 
diagnostic emission lines, for AGNs in the parent, pair, and 
tidal samples (both corrected for fiber incompleteness for the 
latter two), respectively. For these four emission lines, AGNs 
in the pair sample have similar S/N distributions to those of 
the parent AGN sample. On the other hand, AGNs in the 
tidal sample have larger emission-line S/N than do the parent 
AGNs; tidal AGNs have median emission-line S/N which are 
^1.5 times those of parent AGNs. This may be explained if 
AGN pairs with more prominent tidal features tend to exhibit 
stronger emission lines (and thereby have higher S/N detec- 
tions). We increase the S/N threshold and redo the analy- 
sis. In the right panel of Figure [12] we show how the frac- 
tion of AGN pairs (and those with tidal features) varies with 
increasing S/N threshold. The fraction of AGN pairs drops 
from - 3.6% at S/N> 3 to - 1.4% at S/N> 30. Because 
the pair and parent samples have similar S/N distributions, 
the pair fraction would scale ocN^/N (xN if AGNs were un- 
correlated (i.e., randomly paired), where denotes the total 



number of parent AGNs. The obsei-ved decay of pair frac- 
tion with decreasing is less steep than the uncoiTelated case 
would predict. The fraction of AGN pairs with detected tidal 
features roughly stays constant with increasing S/N threshold, 
although we caution that it should be treated as a lower limit 
and the incompleteness due to the inspection of tidal features 
is likely higher at lower S/N, if AGN pairs with more promi- 
nent tidal features tend to exhibit stronger emission lines. At 
S/N> 15, > 80% of all AGN pairs show tidal features (Figure 
[13. 

4. DISCUSSION 
4.1. Implications for AGN Pairs with Smaller Separations 
The fraction of AGN pairs at 0.02 < z < 0.16 with 5 h'^^ 
kpc < rp < 100 /i^Q kpc is '--^3.6% corrected for SDSS spectro- 
scopic incompleteness. The pair fraction decreases to ^1.0% 
if we restrict ourselves to those pairs that show clear tidal fea- 
tures in SDSS images. If we assume that (1) the fraction 
of observed pairs within a certain range of projected sepa- 
rations scales linearly with the time r that a merger spends 
in that range, and (2) the probability that the two AGNs in 
a merger are simultaneously active does not strongly depend 
on the merging phase over the range of separations we are 
considering (i.e., a few kpc to tens-of-kpc scales), then the 
fraction of AGN pairs on tens-of-kpc scales can be related to 
the fraction on kpc scales as /kpc w /tens-of-kpc Tkpc/rtens-of-kpc- 
The timescales r are proportional to the length scales, so /kpc 
will be ^0. 1 % or ^0.4% using the tidal and full pair samples, 
respectively. 

Several studie s (iLiu et all 1201 Obt ISmith et al.l 120101; 
'Wang et al.l2'009h have found that - 1% of SDSS AGNs show 
double-peaked profiles in [Oiii] AA4959,5007 emission lines, 
in which the two velocity peaks are blueshifted and redshifted 
from the systemic v elocity dete r mined from stellar absorption 
fines. As shown in lShen et aH (1201 Ih . at least --50% of the 
double-peaked narrow emission lines are due to narrow-line 
regio n kinematics such as bic onical outflows or rotating disks 
(e^ lAxonetanil998l: IVeilleuxet all 1200 It ICrenshaw et aP 
l2010h . while -10% or more reflect the orbital motion of a 
merging pair of AGNs (e.g. lZhou et al.l2004tlComerford et al.l 
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I2009ah . In the merging AGN pair scenario for those double- 
peaked AGNs, the projected angular separation of the two nu- 
clei has to be smaller than 3" for both AGNs to be covered 
by a single SDSS fiber This corresponds to r,, < 8 h^g kpc 
for the typical redshift of z 0.15 in the ILiu et alJ (1201 Obi) 
sample. The fraction of AGNs in kpc-scale pairs we have in- 
ferred (/kpc ^ 0.1%-0.4%) suggests that only ~ 10%^0% of 
AGNs with double-peaked narrow emission lines are merg- 
ing systems. However, this should be treated as an upper 
limit, considering there are kpc-scale AGN pairs with LOS 
velocity offsets < 150 km s"' and therefore will not be iden- 
tified as SDSS double-peaked narrow-line AGNs. We caution 
that the estimate is highly uncertain since Tkpc/rtens-of-kpc de- 
pends on merger parameters and the host-galaxy properties of 
AGN pairs. In addition, we use the iKauffman n et aL ( 2003a) 
criterio n for diagnostic line ratios, w hereas the Kewlev et al. 
(1200 lb criterion was adopted in the iLiu et al.l (|2010b) sam- 
ple, whose S/N cut for spectral measurem ents was also more 
stringent. Nevertheless, iLiu et al.l (l2010a ) used follow-up ob- 
servations to confirm four kp c-scale bin ary AGNs out of 43 
dou ble-peaked objects of the iLiu et al.l (2010b) sample (see 
also lShen et al.i r201 1). in broad agreement with our estimate 
of /kpc based on AGN pairs with broader separations. 

4.2. Comparison with Binary Quasars 

The fraction of quasars in pairs (with bolometric luminosi- 
ties of > lO'*^ erg s~') with separ ations of tens to hundreds of 
kpc is <0.1 % at 1 < z < 5 (e . g..lHennawi et ani2006l 120101: 
iMyerset all [20081: IShen et al .1120101) . This is ~8-30 times 
smaller than what we observe for the less luminous AGN pairs 
at z ^ 0.08. If we assume that (1) all luminous high-redshift 
quasars are triggered in galaxy-galaxy mergers, and (2) the 
two quasars in a pair shine at random, uncorrected times, 
among all quasars the observed fraction of quasar pairs with 
separations of tens to hundreds of kpc will be roughly the ra- 
tio of the quasar lifetime Tquasm and the time for two galaxies 
with this separation to merge Tmeige- Both Tquasm- and Tmerge 
span orders of magnitude in their parameter spaces. We spec- 
ulate that the first of the two assumptions is likely valid for 
high-redshift luminous quasars, as mergers are arguably the 
most efficient mechanism to give rise to such high-mass ac- 
cretion rates. The second assumption is valid, provided that 
the separations of observed quasar pairs (only a few of which 
are tens-of-kpc separations) are too large for galaxy-galaxy 
tidal interactions to be effective (the counterexample of an 
interacting quasar pair reported by Green et al. 2010 has a 
separation of 2 1 kpc which is at the lower bound of known 
quasar pairs). However, these two assumptions do not nec- 
essarily apply to the less luminous AGN pairs; their intrinsic 
accretion luminosities are roughly three orders of magnitude 
lower, and the typical projected separations in our sample are 
smaller by a factor of ^10 than the high-redshift luminous 
quasar pairs. Our AGN pair sample (low-luminosity, low- 
redshift) and the (luminous, high-redshift) binary quasar sam- 
ples probe populations at very different redshift and AGN- 
luminosity regimes, such that a useful direct comparison with 
the binary quasar samples cannot be made within our sam- 
ple. In particular, the statistics of our sample is dominated by 
Seyfert galaxies (and mostly narrow-line Seyferts; Table |2]l. 

4.3. What Fraction of Moderate-luminosity AGNs Are 
Triggered in Galaxy Interactions? 

As above, if we assume that all moderate-luminosity AGNs 
were triggered in galaxy interactions and the two compo- 



nents in a pair shone at uncorrelated times, the fraction of 
AGNs in pairs with 5 /i^g kpc < < 100 h^^ kpc would 
roughly be the ratio between the typical lifetime of moderate- 
luminosity AGNs tagn and the dynamical timescale of galaxy 
mergers with 5 /i^o kpc < < 100 h^li kpc rs-iookpc, which 
is ^ 10% assuming tagn ^ 10** yr and rs-iookpc ~ 10"^ yr, 
larger than the ^ 1 %-4% value that we have found. This im- 
plies that only - 10%^0%x(10'*yr/rAGN)(T5-iookpc/10''yr) 
of moderate-luminosity AGNs can be triggered in galaxy 
interactions. In addition, we will show in Paper 11 that 
the strengths of AGN activity are correlated between the 
two components in interacting pairs. This correlation sug- 
gests that the two SMBHs in a galaxy pair may tend to be 
activated simultaneously, making the expected fraction of 
AGN pairs even larger than ^ TAON/Ts-iookpc- Thus, the 
observed fraction of AGN pairs suggests that fewer than 
- 10%-40%x(10'*yr/TAGN)(r5-iookpc/10Vr) of moderate- 
luminosity AGNs are triggered in galaxy interactions. 

5. SUMMARY 

We have selected a sample of 1286 AGN pairs (or multi- 
ples) with LOS velocity offsets < 600 km s"' and projected 
separations < 100 h^Q kpc from 138,070 optical AGNs in 
SDSS DR7. Two hundred fifty-six pairs of these show un- 
ambiguous morphological tidal features in their SDSS images 
indicative of ongoing interactions. After correction for spec- 
troscopic incompleteness, the fraction of AGN pairs with 5 
/;-', kpc < rp < 100 /z^jj kpc and Av < 600 km s"' is -3.6% 
among the parent AGNs at 0.02 < z < 0.16. The fraction of 
AGN pairs that show tidal features is —1.0%. 

The current study is a part of our continuing effort to sys- 
tematically identify and characterize the populations of merg- 
ing SMBHs at various stages. The sample that we have pre- 
sented increases the number of known AGN pairs on — 5-100 
hjli kpc scales by more than an order of magnitude, thanks to 
the statistical power of SDSS. It constitutes the starting point 
of a statistical analysis of their properties. In Paper 11, we 
will examine the effects of tidal interactions on AGN pairs 
by quantifying their recent star formation and BH accretion 
activity as a function of pair separation, calibrated against a 
control sample of AGNs matched in both redshift and stellar 
mass distribution, and we will examine correlations between 
the interacting components. 
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APPENDIX 

EXAMPLES OF AGN PAIRS WITH TIDAL FEATURES 

In Figure[T3] we show SDSS images and spectra for three examples of AGN pairs with tidal features. Note that for AGN pairs 
with small (kpc-scale) separat ions (e.g., the top object shown in Figure [TsTl. we cannot rule out the possibility that one AGN 
is ionizing both galaxies (e.g.. iLiu et al.ll2010ah . High-resolution X-ray and/or radio observations can help better constrain the 
double or single AGN nature of these sources. 
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Fig. 13. — Three examples of AGN pairs with tidal features. Left column: SDSS gr('-coIor composite images. North is up and east is to the left. FOV is 
50" X 50" for the top and bottom objects, and 100" X 100" for the middle one. Right column: SDSS fiber spectra. Labeled on each panel are spectroscopic 
plate number, fiber ID, and modified Julian date. Top: 0929-570-52581 is the NW component, and 0930-285-52618 is the SE component. The two AGNs are 
separated by 2. "4 (5.6 h^^ kpc). Middle: the northem component (940-633-52670) is a broad-line, and the southern component (941-255-52709) is a narrow-line 
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Bottom: 1236-54-52751 is the western component which has double-peaked naiTow-emission lines, and 
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